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Considering the possibility to build an e + e~ collider at the energies around the Z°-boson 
resonance with a planned luminosity so high as C oc 10 34 ~ 10 36 cm~ 2 s~ 1 (super Z-factory), we 
make a detailed discussion on the (c&)-quarkonium production through e + + e~ — > (cb)[n\ + b + c 
within the framework of non-relativistic QCD. Here [n] stands for the Fock-states |(cfc)i[ 1 So]}, 
\(f£) a fSo]g), \(db)i?Si]), |(cfe) 8 [ 3 5i] 5 ), KtAJipft]) and Kcb^Pj]) (with J = (1, 2,3)) respectively. 
To simplify the hard-scattering amplitude as much as possible and to derive analytic expressions 
for the purpose of future events simulation, we adopt the "improved trace technology" to do our 
calculation, which deals with the hard scattering amplitude directly at the amplitude level other 
than the conventional way at the squared-amplitude level. Total cross-section uncertainties caused 
by the quark masses are predicted by taking m c = 1.50 ± 0.30 GeV and mj = 4.90 ± 0.40 GeV. If 
all higher (c6)-quarkonium states decay to the ground state B c (|(cfe)i[ So])) with 100% efficiency, 
we obtain o" e + +e -^ Bc+b+s = 5.190^2 419 pb, which shows that about 10 5 ~ 10 7 B c events per 
operation year can be accumulated in the super ^-factory. If taking the collider energy runs 
slightly off the Z°-peak, i.e. y/S = (1.00 ± 0.05)mz, the total cross-section shall be lowered by 
about one-order from its peak value. Such a super Z-factory shall provide another useful platform 
to study the properties of B c meson, or even the properties of its excited P-wave states, in addition 
to its production at the hadronic colliders Tevatron and LHC. 

PACS numbers: 12.38.Bx, 12.39.Jh, 14.40.Pq 



I. INTRODUCTION 



Recently, a Chinese group has raised the proposal of 
constructing a high luminosity e + e~ collider in China, 
the so-called super Z- factory which is similar to the 
Gigaz program suggested by the Internal Linear Collider 
Collaboration but with a even higher luminosity. As 
for this newly purposed super Z- factory, a e + e~ collider 
shall run at the energies around the mass of Z°-boson 
resonance and with a high luminosity up to C cx 10 34 ~ 
10 36 cto _2 s _1 . Many interesting topics can be studied in 
this super Z-factory: the properties of Z°-boson and the 
precision tests of the electro- weak sector of the standard 
model, the properties of the heavy mesons or baryons, r 
physics, the physics beyond the standard model and etc. 

The B c meson has aroused people's great interests 
since its discovery by CDF collaboration H, Q . The B c 
meson is a doubly heavy quark- ant iquark bound state 
and carries flavors explicitly. Systematic studies of its 
production properties at the hadronic colliders Tevatron 
and LHC have been done in the literature, e.g. Refs. 
[l2j . In particular, a computer program BCVEGPY for 
the hadronic production of the B c meson has been com- 
pleted and developed in Refs.[l3l. Il4|. which can be con- 
veniently implemented into PYTHIA [l5] for simulating 
B c events. Very recently, a first observation of the decay 
channel B+ -> J/*7r+7r-7r+ with the help of BCVEGPY 
by LHCb group has been finished At the hadronic 
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colliders, there is much pollution from the hadronic back- 
ground and many produced B c events shall be cut off by 
the trigging condition [UIHI- So, some alternative mea- 
surements would be helpful for a comprehensive study. 

Comparing to the hadronic colliders as the Tevatron 
and LHC, a e + e~ collider has thr ee g eneral features as 
cleanliness, democracy and holism |19| . which are helpful 
and have some advantages to perform more precise mea- 
surements for certain processes. To seek the B c events 
at a e + e~ collider has already been tried by the LEP- 
I experiment, however no B c events have been found 
there. If the luminosity of a e + e~ collider is raised up 
to C cx 10 34_36 cr7i _2 s _1 and its colliding energy is con- 
fined to be around the Z°-peak, it will raise the produc- 
tion rate up to several orders. This is reasonable, since 
as estimated in Ref.[3(, more than 10 9 ~ 11 Z°-boson reso- 
nance can be produced for Gigaz with C cx 10 34 cm~ 2 s _1 . 
Moreover, Refs.(2fjl-|25| have shown that sizable B c events 
can be produce through the Z°-boson decays. Then, as 
a rough combination, one can estimate that such a super 
Z-factory will open new opportunities for studying the 
S c -meson properties. And to be useful reference, it is 
interesting to present a more sound estimation. 

Under the non-relativistic QCD (NRQCD) framework 
[HI , a doubly heavy meson is considered as an expansion 
of Fock states, and the probability of finding each Fock 
state in the expansion is proportional to a definite power 
of v, where v is the relative velocity and v 2 is around 
0.1 ~ 0.3 for heavy quarkonium. We have studied the 
production of (cc)- and (bb)- quarkonium with a photon 
or via initial state radiation at the Z°-peak in the su- 
per Z-factory [27j ■ It is interesting to see whether there 
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has enough B c events produced at such a high luminosity 
e + e~ collider, which can be produced through the process 
e + _+e~ —> 7*/Z° — > (c6)[n]+6+c, where [n]_ stands for the 
(c6)-quarkonium Fock states. Excited (c6)-quarkonium 
states, such as the color-singlet P-wave states and the 
color-octet S'-wave states, may directly or indirectly de- 
cay to the ground state with almost 100% possibility 
via electromagnetic or hadronic interactions. Then, as a 
full estimation of the P c -meson production, we shall dis- 
cuss the following (cfe)-quarkonium states simultaneously; 
i.e. \(cb)^S ]), \(cb) 8 [ 1 S Q }g), |(c&)ipSi]>, \(cb) 8 [ 3 S 1 ]g), 
Kc^xpFr]) and |(c6)i[ 3 Pj]) (with J = (1,2,3)). Here 
the thickened subscript of (cb)-quarkonium denotes the 
color index, 1 for color-singlet and 8 for color-octet; 
the relevant angular momentum quantum numbers are 
shown in the parentheses accordingly. 

Because of the emergence of massive-fermion lines, the 
analytical expression for the squared amplitude becomes 
too complex and lengthy for more (massive) particles in 
the final states and for more higher-order Fock states to 
be generated. For example, to derive the amplitudes for 
the P-wave states, one also needs to get the derivative of 
the amplitudes over the relative momentum of the con- 
stitute quarks. It has been found that to do the numeri- 
cal calculation using the conventional squared-amplitudc 
technology becomes time-consuming for these complex 
processes, since the cross-terms increases with the incre- 
ment of Feynman diagrams, \M\ 2 = J2ij MiM*, where i 
stands for the number of Feynman diagrams for the pro- 
cess. One important way to solve this is to deal with the 
process directly at the amplitude level. For the purpose, 
the "improved trace techn ology " is suggested and devel- 
oped in the literature [H, l21r423l [28j j . After generating 
proper phase-space points, one first calculate the numer- 
ical value for the amplitudes, and then sum these values 
algebraically and square it to get the squared amplitude, 
\M.\ 2 = | J2;-Mi\ 2 ; through such way, numerical simula- 
tion efficiency can be greatly improved in comparison to 
the usual squared amplitude technology Q- Moreover, un- 
der the approach, many simplifications can be done at the 
amplitude level due to the fermion-line symmetries and 
the specific properties of each heavy-quarkonium Fock 
states, then, we can even written down the analytic ex- 
pressions for the amplitude. 

The "Improved trace technology" has been adopted to 
derive a simplified and compact analytical expressions for 
the on-shell Z°-boson decay process, Z° — > (cb) [n] + b + c 
[2TM23I ] . At the present, we will adopt the same method 
for dealing with the (c&)-quarkonium production through 



1 One can also generate B c events through PYTHIA by e+e — an- 
nihilation. However, in PYTHIA all kinds of 6-hadron ( B s , B c , 
A;, and etc. ) events will be produced according to the proper 
fragmentation possibilities, among which the fragmentation pos- 
sibility for the doubly heavy B c meson is quite small [5f|, thus 
using PYTHIA to generate B c events is not 'economical'; i.e., 
when simulating the B c events, too much time will be wasted by 
generating enormous unwanted 6-hadron events. 




FIG. 1. Feynman diagrams for the process e + e — > 
J*/Z° — > (cb)[n\ + b + c, where [n] stands for the Fock 
states: |(c6)i[ f S ]), \(dS)^8a]g), |(c6)i[ 3 Si]>, \(db)sfSi]g), 
Kc&hpPi]} and \(cb)i[ 3 Pj]) (with J= (1,2,3)) respectively. 

the full production process, e + +e~ — s- 7*/Z° — > (c6)[n] + 
b + c. The intermediate Z°-boson may or may not be on- 
shell, and if it is off-shell, some new Lorentz structures 
must be introduced to ensure a gauge- invariant result. 
In a sense, this paper provides an extension of the "im- 
proved trace technology" on how to deal with the off-shell 
particles. 

The remaining parts of the paper are organized as fol- 
lows. In Sec. II, we present the formulation for dealing 
with the process of e + + e~ ->• Z a /j* [cb)[n] + b + c, 
where improved trace technology are adopted to simplify 
the hard scattering amplitude. In Sec. Ill, we give the 
numerical results. Sec. IV is reserved for a summary. All 
necessary expressions are put in the Appendices. 

II. FORMULATION AND CALCULATION 
TECHNIQUES 

As shown in Fig.fl]), there are four Feynman dia- 
grams for the process e + (p2) + e^(pi) — > Z°/j* — > 
(cb)[n](q 3 ) + b(q 2 ) + c{q\), where [n] stands for the 
(c6)-quarkonium Fock-states: J(c6)i[ 1 S'o]), Kcb)^ 1 So]g) , 
Kcfijxpft]), |(c6)8[ 3 Si]3>, \(cb) 1 [ 1 P 1 )) and Kcb^Pj}) 
(with J = (1,2,3)) respectively. Here as a comparison, 
the production channel via a virtual photon, which is 
small around the Z° peak, is also presented. 

According to NRQCD factorization formulae [l^, the 
cross-section of the process can be written in the follow- 
ing factorization form: 

da = ^da{e+ + e~ -> (cb)[n] + b + c)(0 H (n)), (1) 

n 

where the matrix element (O (n)) is proportional to the 
inclusive transition probability of the perturbative state 
(c6)[n] into the bound states. da(e + +e~ — > (cb)[n]+b+c) 
stands for the short-distance cross-section; i.e. 

d&(e + + e~ -> (cb)[n] + b + c) 

= — \ ==g|M| 2 d$ 3 , (2) 
4\/ (Pi ■ P2) - m\m^ ^— ' 
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where ^2 means we need to average over the spin states 
of initial particles and to sum over the color and spin of 
all final particles, and the three-particle phase space 
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The phase space can be generated by the routine RAM- 
BOS [30], and its integration can be carried out numeri- 
cally with the VEGAS program (3l| . 

The scattering amplitude for the process can be writ- 
ten as, 



where 



and 



iM =C x L^D^H^,, 



^ r ,=«rfe)r^(pi) 



(3) 



(4) 



(5) 



The spin- indices r, r' , s and s' stand for the spin pro- 
jections of the initial electron and positron, and the 
outgoing quark and antiquark respectively. The over- 
all parameter C = C s or C is for the production of 
the color-singlet or the color-octet quarkonium state re- 
spectively. For the quarkonium production through Z°- 
boson propagator: C s = ^= JL gk_ anc { Q = 

■Tj- x {y/2T a T b T a ) lJ (V2T b stands for the color 



sin 2 S„ cos 2 S„ A \ vaJ j; - 1 

of the color-octet (c6)-quarkonium state); the vertex 
pp _ _ 4 s in 2 — 7 5 ) and the propagator = 
k^-m%+im z v, (~9iw + k^kv/k 2 ) , where T z stands for the 
total decay width of Z° boson. If the e + e~ collider runs 
at the energies around the Z° mass, it will constitute 
an energy scan. The shape of the cross-section varia- 
tion around the Z° peak is described by a Breit-Wigner 
ansatz with an energy- dependent total width. While, for 
the quarkonium production through photon propagator: 
C s = ^e 2 5 2 % and C = e 2 g 2 (V2T a T b T a )if, the vertex 
pp = y an( j the propagator D^ u = j^-g^v The gamma 
structures A£ with n — 1, • • • , 4 can be read off from the 
four Feynman diagrams in Fig.([T]), whose explicit expres- 
sions can be directly found in Refs.|22l [23j| q 

Ly r , and H v ss , contain explicit strings of Dirac 7- 
matrices sandwiched by two fermion spinors, which can 
be done under the conventional trace technology. How- 
ever it is hard to write down the squared amplitudes for 



these massive processes, especially for the P-wave cases, 
which is too complex and lengthy. Here, to simplify the 
analytic results and to make it more compact, we adopt 
the "improved trace technology" [12i l2lW23l | . It is noted 
that under suitable transformations, this approach is also 
useful for deriving the polarized cross-sections, which is 
helpful for a polarized e + e _ collider. 

As an explanation for the "improved trace technology" , 
we first arrange the amplitude H" s , into four orthogonal 
sub-amplitudes H" according to the spins of the outgo- 
ing quark q with spin s and antiquark q with spin s', 
then transform these sub-amplitudes into a trace form 
by properly dealing with the massive spinors with the 
help of an arbitrary light-like momentum fco and an ar- 
bitrary space-like momentum ki, which satisfies k\ = — 1 
and fco • fci = 0. The final results are independent of 
fco and fci, which provides a way to check the Tightness 
of the derived results. And one can choose them to be 
those that will maximumly simply the amplitude. And 
then we do the trace of the Dirac 7-matrix strings at 
the amplitude level, which shall result in explicit series 
over some independent Lorentz-structures. The ampli- 
tude L^ r , can be done in a similar way. Detailed process 
of the approach can be found in Refs.[22|,[23j], so to short 
the paper, we only list all the differences for the present 
processes. However, because at the present Z° may or 
may not be on-shell, one should use these on-shell ex- 
pressions carefully with suitable transformations. 

As for H" s , , it can be rearranged as the following four 
spin combinations 



(6) 
(7) 
(8) 













H v 2 = 


H's* 


Hiss' 






#3 = 




s> - Hiss' 






Ht = 




S ' + H1 ss , 




V2 



(9) 



The production through the photon propagator is simi- 
lar to the case of the Z°-boson propagator, which can be 
directly read from its vector part. So we shall only con- 
centrate our attention on the case of the Z°-boson prop- 
agator 0. iff 2 3 4 can be expressed over the independent 
Lorentz structures. For example, for the production of 
I (c&)i[ #0]) -quarkonium, we have 



So1 1 m z 2 m z 3 m z 
+A[ Sa] -^-e(iy,k,q2,q 3 ), 



(10) 



2 As for the case of Z° propagator, from the correspondence be- 
tween A n and An, one only needs to replace ^(fc) in the ampli- 
tudes of Refs. [23,1231 by 7". And the case of photon propagator 
can be easily derived by changing the corresponding vertex. 



3 A generator BEEC for dealing with the heavy quarkonium pro- 
duction at a e + e - collider within the PYTHIA environment [lf| 
is in preparation, where both cases with the photon and the Z°- 
boson propagators are included |3'J . 
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where e(i/, k, 92,93) = e v ^ pa kpq^pq^a- Comparing the 
Lorentz expansions of H" 1S( ^ with those of Mp.s ] m 
Refs.[22|, we obtain the same expansion coefficients 
for the corresponding Lorentz structures. More explic- 
itly, by multiplying both side of Eq.lflO]) with the on- 
shell Z°-boson polarization vector e„(fc), which satisfies 
k • e(k) = 0, we can obtain exactly the same matrix el- 
ements i.e. e 1/ (fc)i?n So j|fc. e= o = M[i5 j. However, when 

the Z°-boson is not on-shell, (A^ So] k v /mz) must be 
kept in HKg,. And to be a consistent analysis for the 
present case, we need to take these new terms into consid- 
eration. For convenience, we present all the new Lorentz 
structures for all the mentioned Fock states together with 
their non-zero coefficients in Appendix B. 

As for the leptonic vector L^ r ,, it is simpler and can 
be directly expressed as 



M M A* 

rr ' 1 m z 2 m z 3 m z 



(ii) 



where all the non-zero coefficients C^ r are presented in 
Appendix C. 

To be a cross check, we also use the squared amplitude 
approach to deal with the process, which however can be 
rewritten as the following more convenient form 



1 



C 2 x (D a/3 L a P) ■ (D pa H i 



pa J 



The leptonic tensor 

L a0 = _ 4g «P pi . p2 + 4p /p 2 a + 4 Pl >/, 



(12) 



(13) 



where small mass terms involving m e are neglected. The 
result for the phase-space integration over {D pa H pa ) is 
equal to the decay width of Z° — > (cb) [n] + b + c times a 
factor l/[(fc 2 — m 2 z ) 2 + m%F%\ . Numerically, we obtain 
the same results for these two approaches. 



III. NUMERICAL RESULTS 

The basic input parameters are chosen as the follow- 
ing values: [33]: m z = 91.1876 GeV, T z = 2.4952 GeV, 
sinX = 0.23119 and a s (m z ) = 0.1176. We adopt 
the leading-order a s running, and take the normaliza- 
tion scale to be 2m c , which leads to a s (2m c ) = 0.212. 
The color-singlet matrix elements can be directly related 
to the wave functions at the origin for the 5-wave states 
or with the first derivative of the wave functions at the 
origin for the P-wave states. The S'-wave color-octet 
ones are suppressed to the corresponding ones by a factor 
A 2 s (v) with A s (v) 



im,l26j,l3J]. In the leading order 
approximation, the color-octet or color-singlet matrix- 
element appears in the amplitude as a linear factor. The 
uncertainties in the production from the matrix elements 




E(GeV) 

FIG. 2. Total cross section a(e + + e~ -> -y*/Z° -> (cb)[n] + 
b + c) versus the e + e~ collision energy E, where contributions 
from the mentioned Fock states |(c6)i,8[ 1 So]}, j (c&^spSi]}, 
|(c6)i[ 1 Pi]} and |(c6)i[ 3 Pj]} have been summed up. Here the 
dashed, the dotted and the solid lines are for the production 
channel with 7*-propagator, Z°-propagator and the sum of 
these two channels, respectively. 



can be figured out straightforwardly, so throughout the 
paper, we will fix the values of the radial wave function 
at the origin and the first derivative of the radial wave 
function at the origin to be |Ps(0)| 2 = 1.642 GeV 3 and 
|P' P (0)| 2 = 0.201 GeV 5 [H which are derived under the 
Buchmuller-Tye-potential model [36[ . 

To calculate the hard-scattering amplitude for the 
short-distance cross-section da, the equations: q 2 = m 2 , 
q? = ml and q 2 = M 2 (M is the (c6)-quarkonium mass, 
93 = 9c + qi with 9 C and qi are momenta of the constitute 
quarks) must be satisfied simultaneously, and only when 
these equations are satisfied can the gauge invariance of 
the amplitude be guaranteed. In this paper, we proceed 
the numerical calculation by taking m c = 1.50 GeV and 
nib = 4.90 GeV as their center values @- The mass of 
the quarkonium bound-state is the sum of the two heavy 
quark masses; i.e. M = nit + m c = 6.40 GeV. We will 
discuss the theoretical uncertainties caused by the c- and 
b- quark masses by varying m c = 1.50 ± 0.30 GeV and 
m b = 4.90 ± 0.40 GeV. 

Firstly, we calculate the total cross-section for the 



"e collider energy at the Z° peak, i.e. E = m z . We 



obta 



^(cbKPSo]) = 1-800 pb, 
a \(cb)i[ 3 Si}) = 2.575 pb, 
CKcfcHppj) = 0.190 pb, 
"Kc&HpPb]) = 0.115 pb, 
°l(c5)ippi]) = 0- 232 P b > 
fKc&HPft]) = 0.256 pb, 



(14) 
(15) 
(16) 
(17) 
(18) 
(19) 



4 There values are consistent with the so called IS masses |37|| . i.e. 
half of the J/^-mass for charm quark and half of the T(1S) for 
bottom quark. And they are also consistent with the bound-state 
parameters under the potential-model calculations [35l l3SI| . 
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FIG. 3. The differential pt- and y- distributions of the 
(c5)-quarkonium, where the solid, the dashed, the dotted, 
the dash-dot, the triangle and the diamond lines are for 
KcbhtSo}), \(cb) 1 [ 1 P 1 ]), !( C &)i[ 3 A]>, |(cb) 1 [ 3 P 2 ]), j(c6)i[ 3 Si]> 
and |(c&)i[ 3 Po]} respectively. 





95% m z 


97% m z 


m z 


103% m z 


105% m z 


0tot (pb) 


0.341 


0.860 


5.190 


0.930 


0.384 



TABLE I. Total cross-section of the process e + + e~ — > 
7* /Z° — > (cb)[n]+b+c under different collision energies, where 
contributions from the mentioned Fock states |(c6)i,b[ So]), 
\{d))iA 3 Si]), \{cb) 1 [ 1 P 1 ]) and \(cb) 1 [ 3 P J ]) have been summed 
up. Here m c = 1.50 GeV, m b = 4.90 GeV and v 2 = 0.2. 



<r\W.?s )g) = 0-225 x« 4 pb, (20) 

^|(o5).ps 1 ] 9 ) = 0.322 xv 4 pb. (21) 

It is found that the total cross-sections of all the P-wave 
states can be added up to ~ 20% of that of the ground 
state i? c (|(c6)i[ 1 S'o]). So the P-wave (c6)-quarkonium 
states should be taken into consideration, especially for 
the future high luminosity e + e~ colliders. Two color- 
octet iS-wave states shall only provide ~ 1% contribution 
to the ground state. 

As a reference, we present the total cross section 
cr(e+ + er -> j*/Z° -> (cb)[n] + b + c) versus the e + e" 
collision energy E in Fig. ([2]), where contributions from 
the mentioned Fock states {(cfyx^So]), |(c6)i i8 [ 3 <S'i]), 
|(c5)i[ 1 Pi]) and |(c&)i[ 3 Pj]) have been summed up. We 
take v 2 = 0.2 for the color-octet states. For small e + e~ 
collision energy (E < 60 GeV), we have 

Cr e++e--i-7*->(cfo)[n] > <7 e++e- -^Z n ^(cb) [n] - 



m c (GeV) 


1.20 


1.35 


1.50 


1.65 


1.80 


(7 l(=5)i[ 1 Sol)(P b ) 


3.564 


2.485 


1.800 


1.477 


1.031 


<7 l(c5)i[ 3 Si]>(P rj ) 


5.447 


3.668 


2.575 


2.074 


1.397 


cr l(^)if 1 Pil)(P b ) 


0.571 


0.319 


0.190 


0.139 


0.079 


cr l(=5)ipP l>(P b ) 


0.281 


0.175 


0.115 


0.089 


0.056 


cr l(=5)i[ 3 f > i])(P b ) 


0.658 


0.379 


0.232 


0.172 


0.101 


cr l( C 5) 1 pP 2 l)(pb) 


0.845 


0.450 


0.256 


0.154 


0.097 



TABLE II. Uncertainties for the total cross-section of the pro- 
cess e + + eT — > (cb)]n\ + b + c with varying m c , where mj is 
fixed to be 4.90 GeV. 



m b (GeV) 


4.50 


4.70 


4.90 


5.10 


5.30 


CT l(ci>) 1 [ 1 Snl)(P b ) 


1.816 


1.808 


1.800 


1.793 


1.786 




2.526 


2.551 


2.575 


2.599 


2.623 


(T l( C S) 1 [ip 1 l)(P b ) 


0.195 


0.192 


0.190 


0.188 


0.186 


(T l( C S) 1 f3p l>(pb) 


0.126 


0.120 


0.115 


0.110 


0.105 


t7 l(=S)ipJ=il>(P b ) 


0.243 


0.237 


0.232 


0.227 


0.223 


cr l( C 5) 1 pP 2 l>(P b ) 


0.252 


0.254 


0.256 


0.258 


0.260 



TABLE III. Uncertainties for the total cross-section of the 
process e + + eT — ► (c6)[n] + b + c with varying m b , where m c 
is fixed to be 1.50 GeV. 



One may observe a small peak around 20 GeV, which is 
caused by the production channel through 7* -propagator 
and has already been observed by Ref.[21|. When the 
collider energy runs from ~ 10 GeV to 100 GeV, it is 
found that the total cross-section at the energy around 
Z°-peak raises up to four orders than that of 10 GeV due 
to the Z°-resonance effects. Moreover, as shown by Table 
HI if taking the collision energy E = (1 ± 3%)mz, the 
total cross-section shall be lowered to ~ (17%) of its peak 
value; and if taking the collision energy E = (l±5%)mz, 
the total cross-section shall be lowered to ~ 7% of its 
peak value. Moreover, we present the (cb)-quarkonium 
transverse momentum {jpt)- and rapidity (y)- differential 
distributions in Fig. (J3J) . 

Secondly, we discuss the theoretical uncertainties 
caused by the c and b quarks by varying m c — 1.50±0.30 
GeV and m b = 4.90 ± 0.40 GeV. There are other uncer- 
tainty sources include the matrix elements, the renormal- 
ization scale and etc. There are overall factors, their 
uncertainties can be easily estimated when their values 
are known, so we shall not discuss them there. The cross- 
section for the (c6)-quarkonium production through e + e~~ 
annihilation with varying m c and m b are presented in 
Tabs. . It can be found that the total cross-section 
is more sensitive to m c than that of m b . And the to- 
tal cross-section of the P-wave states are more sensitive 
to quark masses than the 5-wave states. By adding the 
uncertainties in quadrature, we obtain 

*|(c5)xrs ]> = l-800tS;?& pb, (22) 
(7 \(cb) 1 [ 3 s 1 ]) = 2.575j:^7g pb, ( 23 ) 
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Pt(GeV) 




FIG. 4. Uncertainties of pt- and y- differential distri- 
butions for the (cfo)-quarkonium, where contributions from 
the Fock states Kc&KspSo]}, |(c6)i, 8 [ 3 A]>, Kc^i^Pi]) and 
|(cfe)i[ 3 Pj]} have been summed up. 





'-Pi]) 


= 0.190±g;ffi 


pb, 




(24) 




3 Po]) 




pb, 




(25) 




3 Pi]) 


n 9 no+0.426 
— U.zoz_ 131 


pb, 




(26) 


CT IM>)i 


3 P2]) 


= 0.266l8;Sg 


pb, 




(27) 




s )g) 


n nnr+0.221 

— u.zzo_ 096 


4 

X t> 


pb, 


(28) 




Si]g) 


n qoo+°- 359 
— u.ozz_q -^47 


4 

X V 


pb. 


(29) 



If assuming the higher excited states decay to the ground 
state B c (|(c6i)[ 1 5'o])) with 100% efficiency, then we ob- 
tain the total cross-section for B c , i.e. 



e++e~ 



ci qo+6.222 

: +&+c — °- iJU -2.419 



pb. 



(30) 



Furthermore, to show the corresponding uncertainties 
more clearly, we present the uncertainties of pt- and 
y- differential distributions for the (rf>)-quarkonium pro- 
duction in Fig. ((H), where the Fock states |(c6)i s^So]), 
|(c6)i, 8 [ 3 Si]), Kc^ipPx]) and Kcb^Pj}) states have 
been summed up. The uncertainties are drawn by a 
shaded band, where the center solid line is for m c = 1.5 
GeV and mi, = 4.9 GeV, the upper edge of the band is 
for m c = 1.2 GeV and m& = 5.3 GeV, and the lower edge 
of the band is for m c = 1.8 GeV and nib = 4.5 GeV. 

Experimentally, there is no detector which can cover 
all the kinematics of the events, so only some of the B c 
events can be observed. That is, the B c events with small 




FIG. 5. The (cfo)-quarkonium differential distributions versus 
pt with various y C ut, where contributions from all the con- 
sidered Fock states have been summed up. The solid, the 
dashed, the dotted, the dash-dot and the dash-dot-dot lines 
are for y cnt = 0, 0.5, 1.0, 1.5 and 2.0 respectively. 




FIG. 6. The (cfe)-quarkonium differential distributions versus 
y with various ptcut, where contributions from all the con- 
sidered Fock states have been summed up. The solid, the 
dashed, the doted and the dash-dot lines are for ptcut = 0, 5, 
20, 35 GeV respectively. 



Pt and/or a large rapidity y cannot be detected, so they 
cannot be utilized for experimental studies. As a final 
remark, it would be helpful to know how the differen- 
tial cross-sections change with different pt- or y- cuts. In 
Fig. (JS|) , we present the (c6)-quarkonium ^-distributions 
for all the considered Fock states with a particular rapid- 
ity cut (y cu t), where the solid, the dashed, the dotted, 
the dash-dot and the dash-dot-dot lines are for y cut — 0, 
0.5, 1.0, 1.5 and 2.0 respectively. In Fig.®, we present 
the (c6)-quarkonium ^-distributions for all the consid- 
ered Fock states with a particular transverse momentum 
cut (ptcut), where the solid, the dashed, the doted and 
the dash-dot lines are for p tcu t — 0, 5, 20, 35 GeV respec- 
tively. To analyze the quantitative difference of the differ- 
ential distributions with regard to the parameters ptcut 
and y cu ti we introduce a ratio for the integrated total 
cross-section: 



FL 



(31) 
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TABLE IV. R Ptcut versus the parameters ptcut and the y c 



Ptcut 


GeV 


5 GeV 


20 GeV 


Vent 


0.5 


1.0 


1.5 


0.5 


1.0 


1.5 


0.5 


1.0 


1.5 




0.382 


0.710 


0.897 


0.391 


0.725 


0.913 


0.520 


0.913 


1.00 
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Appendix A: Linear polarization vector and tensor 



where cr ycMt and <To are the total cross-section with and 
without y cut respectively. In the definition of R ptaut , 
contributions from all the considered Fock states have 
been summed up. The values of the ratio P Ptc „ t that 
correspond to different values of ptcut and y cu t are pre- 
sented in Tab.JlVl). It is found that when under the same 



Ptcut, Rp tcut increases with the increment of y C ut, and 
when under the same y cu t, Rp taut also increases with 
Ptcut- This shows that the rapidity distribution of the 
(c6)-quarkonium shrinks with the increment oiptcuu e.g. 
when ptcut = 20 GeV, \y\ < 1.5. 



We set the z-axis to be along the incoming 
beam. The momentum of the (c&)-quarkonium qt^ = 
(q g, \q 3 \ sin6>co s</>, \q 3 \ sinflsin^, \q 3 \ cosfl), where \q 3 \ = 
V {Q3) 2 ~ M 2 > M is the quarkonium mass, 8 and (j> are 
polar and azimuth angles respectively. The expressions 
of the polarization vector can be written as 

e x(q 3 ) = (0, cos£? cos 4>, cos 6* sin0, — sin (9), 
£ y(?3) = (O,-sin0, cos0, 0), 

£z(q 3 ) = —Qq%\,q^8ia9coB(f>,q^sm.e8ia<f>,q^ cos 6). 

As required, the polarization vector satisfies the following 
relation 



IV. SUMMARY 

In the present paper, we have studied the (cb)- 
quarkonium production through e + e~ annihilation in 
leading a s order, but with the u 2 -expansion up to 
v 4 . The production of the (c6)-quarkonium in Fock 
states Kc&JxpSo]), l(c&)i[ 3 Si]), \(cb) 1 [ 1 P 1 }) : {(ob^Pj]), 
\(cb) 8 [ 1 S ]g) and |(c6) 8 [ 3 S' 1 ]g) have been studied. Higher 
Fock states such as the P-wave states will also provide 
sizable contributions to the total cross-section of the B c 
meson. If assuming the higher excited states decay to 
the ground state i3 c (|(c6i)[ 1 S'o])) with 100% efficiency 
via electromagnetic or hadronic interactions, then we will 
obtain a. 



e t ■ q 3 = , e, • ej 



(i,3=x,y,z). (Al) 



e++e-->B c +b+E 



5.190±|;2?| P b. 



According to our present calculation, one would ex- 
pect to accumulate 10 5 ~ 10 7 B c events in one opera- 
tional year for the super Z-factory with a high luminosity 
C oc 10 34 - 10 36 cm~ 2 s _1 . If taking the e + e~ collision 
energy to run slightly off the Z°-peak; i.e. y/S = 0.95m^ 
or 1.05m^, the total production cross-section will be low- 
ered by about one-order from its peak value. However one 
can estimate 10 4 — 10 6 B c events to be produced under 
such condition, which is still sizable. Then, in addition 
to the hadronic colliders Tevatron and LHC, the super Z- 
factory will provide another good platform for studying 
the P c -meson properties. 



The polarization tensor €^(q 3 ) for ( 3 P2)-quarkonium 
state can be constructed in terms of e x (q 3 ), e y (q 3 ) and 
£z{q 3 ) as 



J-( f V f v _|_ fPf"} 

x V y x) s 
fc 3 — V * ~ z y) 1 

71 v x x 



uv _ 



2 



,/^\ c x c x c y c yJ 



V2 



(A2) 



Appendix B: New lorentz structures and their 
coefficients for e + + e~ — > Z° — > (cb)[n] + b + c 

In the present section, we present the new lorentz 
structures and their non-zero coefficients for the process 
e+(p 2 ) + e"(pi) -> Z°/Y -> (cb)[n](q 3 ) + b(q 2 ) + g(«i), 
where n stands for Kcb)^ 1 S }) , \(cb)i^Si]), \(cb)t?Pi]), 
|(c6)i[ 3 Pj])(J=l,2,3) respectively. All these new Lorentz 
structures have no correspondence to the Lorentz struc- 
tures introduced from the process Z°(k) — > (cb)[n](q 3 ) + 
b{q%) + c(qi). Note, the following parameters, r\, r 2 , 
r 3 , Tb, T c and so on, have the same meaning as those 
of Refs.[22|, [23|, and we shall not present them here for 
short. 



1. Spin-singlet S-wave state: |(c6)[ x Sb]) 
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in part by the Fundamental Research Funds for the 
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There is only one new Lorentz structure for |(c5)[ 1 S'o]), 

k v 



Pi = 



m z 



(Bl) 



whose coefficients A*' 2 ' 3 ' 4 are 



A* 



m z 



7 / 2 Li 



(d 4 ri(3ri — r 2 )r 3 + di(3r 3 ri 2 + (— r 2 r 3 + u — 2x)n — r 2 u) 



+(d 2 (r 1 (3r 1 + r 2 ) + u - 2y) + d 3 (ri(3ri + r 2 ) + 3u - 2y))r 3 + d 4 (3r! + r 2 )u), 



4m z 7/2 L: 



(((d 3 + d 4 )T c - d 2 T b )(n{3 ri - r 2 )r 3 + (3n + r 2 )u) + diT 6 (3r 3 r 2 - (r 2 r 3 



-3u + 2ar)n + r 2 u) + 2r 3 (d 2 T b - d 3 T c )y) 

(diT 6 (u(2y - 1) + 2nr 3 (x + 2y - 1)) + d 2 T b (u - 2(-r 2 r 3 y + nr 3 (x + y - 1) 
))) + T c (2di(r ir3 + u){2nr 3 + u) + d 3 (u(2y - 1) + 2nr 3 (x + 2y - 1)))) 
(m z 9/2 N (2d i (2r* + r 2 r\ + 3urj + r£n + u(u - r 2 )) - d 3 (-4yr 2 



+2(r 2 r 2 - (r| + x - l)n + r 2 ?i)r 3 + u - 2uy) + d 2 (« + 2(r? + (-r 2 . +u + x 

+y - l)r\ + r 2 (u + x + 2y - l)n + r\y - u(x + y))) + di{&r 3 r\ 

-2(-3u + 3x + 2y- l)rj - 2r 3 r 2 r 2 + 2r 2 (u - x - 2y + l)n + u(2y - 1)))) 



(B2) 



(B3) 



(B4) 



2. Spin-triplet S-wave state: |(c6)[ 3 Si]) 
There are three new Lorentz structures for |(c&)[ 3 Si]), 



in 



1 1 

-rk-e(q 3 )k u , B 2 = —q 2 ■ e(q 3 )k u ', 



B 3 = —e(k,q 3 ,q 2 ,e(q 3 ))k v , 



ml 



(B5) 



whose coefficients A 1 ^ 2 ' 3 ' 4 (j = 1,2,3) are 



A\=%d Y L imz ,l2 r 1 ^ r 3 T h , (B6) 
A\ = 4L imz y 2 r 3 3 /2 [(d 2 - di)T 6 - (ds + d 4 )T c ], (B7) 



! 2L 1 mz 7 l 2 fii 
A 3 = — — (d 2 - d 3 ), 



(B8) 

A 2 = 2d 1 i 2 m z 7 / 2 r 1 ^ 3 -, (B9) 
A\ = -L 2 m z r/2 y / : F 3 '[(3d 1 -d 2 + d 3 + d 4 )n 

-(d 1 -d 2 +d 3 -d 4 )r 2 ], (BIO) 



A2 = 8L 2 m z 7 / 2 



(d 2 T b + d 3 T c ), 



. 3 > dunz^Npy/^ 
A i = — (nr 3 + u), 



(Bll) 



(B12) 



3 / m z 9/2 N ^ 2 
^2 = — — - — (d 2 - d 3 + 2d 4 (r 3 + u) 

4:L 2 

+d 1 (-l + 2rlr 3 + 2x + 2y) + 
2(d 2 (r 1 r 3 - x - y) + d 3 (r 2 r 3 + y))), 



, _ 2m z 9 / 2 N ^ 



{d\T b + d 3 T c ), 



Ai , 2d imz 9 / 2 N 0v ^T b 2 , 
A\ = {r( - r 2 n + u), 

.4' rn z 9/2 N 0y /r^ 

A 2 = — (d 2 T b {-2x - 2y + 1) + 

diT b (2x + 2y - 1) + T c (d 4 (rf - r\ + r§ 



+2«) + ds(2l/-l))), 



1 4' 



mz^Vo 
2^1^ 



(B13) 
(B14) 

(B15) 



(B16) 



[(<fa - ds)(ri - r 2 ) + (d2 - di)r#17) 



3. Spin-singlet P-wave state: |(c&)[ 1 Pi]} 

There are three new Lorentz structure for |(c6)[ 1 Pi]), 
which are the same as that of (cfe)[ 3 «Si], whose coefficients 

A)' 2 ' 3 '' 4 ' (j = 1,2,3) are 



yy^ 5/2 

A\ = - F J-{2d xl r 3 r x r\ - 2((d 2 i + d 22 - ci 32 )ri(3ri + r 2 )r 3 + (d 2 i + d 22 - 3d 32 )ur 3 



+dn(3r 3 r 2 + ur\ - 2xr Y - r 2 u))r 2 + 4(d 2J + d 22 - d 32 )r 3 yr 2 + di(r 2 - ri)r 3 ), (B18) 
A 2 = — =(m z 5/2 L 1 (-d 1 r 3 r 2 + d 4 r 3 r% + d 2 r\r 2 - d 3 r\r 2 - 8(d 22 + d 31 )nr 3 yr 2 

+4d 11 r 1 (3r 3 r 2 + (u — 2x)r x - r 2 u - r 2 r 3 r{)r 2 + r 1 (-4d 41 r 3 r 1 r 2 
+4(d 4 i(3ri + r 2 ) + (d 22 + 3d 3 i)r 3 )w 2 + r 3 (4(d 22 + d 3 i)rir 2 

-2(-6(d 22 + d 3 i + d 4i )r 2 + di + d 4 )r 2 + (3di + d 4 )r 4 + (d 3 - d 2 )r 3 )))), (B19) 



^ 2 {di(rl - r{)T b + 2r 2 (((d 21 + d 22 )T b + d 32 T c )(3r 3 r 2 + 3t*n + r 2 u - 2r 3 y 

-r 2 r 3 n) - diiT 6 (3r 3 r 2 + 3t*n - 2an + r 2 u - r 2 r 3 n))), (B21) 

^— (2L 2 m z 5 / 2 (rf 1 r 3 T b (r 1 - r 2 f - d 2 r 3 T b (n - r 2 f + d 3 r 3 T c (n - r 2 f + d A r 3 T c {n 

r\r 2 yjri 

-r 2 ) 2 + 4diirir 2 T 6 (3r 3 r 2 + 3uri - 2xri + r 2 u - r 2 r 3 ri) - 4:d 22 rir 2 T b (3r 3 rl + iur 1 
+r 2 u - 2r 3 y - r 2 r 3 r x ) + 4d 4 ir :L r 2 X' c (3r 3 r 2 + 2,ur x + r 2 u - r 2 r 3 ) 

+4d 31 r 1 r 2 T c (3r 3 r 2 + 3un + r 2 u - 2r 3 y - r 2 r 3 ))), (B22) 

mzb ' 2 ^ L \ d 2 -d 3 ), (B23) 
r\r 2 

l -^{m z 7 ' 2 N {d 1 {r 2 - ri )r 3 T b { ri r 3 + u) + 2d 21 r 1 r 2 T b (2r 1 r 3 (x - 1) + u(2x - 1)) 



+2d 22 r 1 r 2 T b (2r 1 r 3 (x - 1) + u(2x - 1)) + 4(d 2i + d 22 )ri(n - r 2 )r 2 r 3 T b y + 4(d 2 i + d 22 )nr 2 
xT b uy + 2dnr 1 r 2 T b (—2yu + u — 2r x r 3 (x + 2y - 1)) 

+2d 32 r 1 r 2 T c (u(2y - 1) + 2r 1 r 3 (.x + 2y - 1)))), (B24) 

— —(m z 7/2 N (d 3 r 3 (r 3 (r 3 - r 2 )r 2 + r 2 + r x {r 2 r 3 - l))T c + d A {r x - r 2 )r 3 {r\ - r\ 

2L 2 rir 2 Jri 



+r 2 3 + 2u)T c + Ad 41 r 1 r 2 (2r 1 r 3 + u)(r{ -r\+r\ + 2u)T c + 4d 3 irir 2 (2rir 3 (a; - 1) - u)T c 
+2d 32 r 1 r 2 r 3 {2r 2 x - r 3 )T c + 2(d 3 (r 1 - r 2 )r 3 + 4rf 3i rir 2 (2rir 3 + u))yT c + 2d 21 r 1 r 2 r 3 T b (r 3 
-2r x x) + di(n - r 2 )r 3 T b {2x + 2y - 1) - M xx r x r 2 T b (-2yu + u- 2 ri r 3 (x + 2y- 1)) 
-d 2 r 3 T b {2r 2 r\ + (2r 2 + 2x + 2y - l)n + r 2 - 2r 2 (a; + y)) 

+4d 22 r 1 r 2 T fc (« - 2(-r 2 r 3 y + ri r 3 (x + y - 1) + u(a; + y))))), (B25) 



(-2d 2 n +di(n -r 2 )+d 3 r 3 ) (B26) 
4rir 2 L 2 

7/ 2 

. Z ° (-4dnr 2 r 3 r 2 rf + 2d n r 2 (-u + 2n(r 3 + (3n + r 2 )(u - x) + 3) 
4rir 2l /r^Li 

+2(w - 2rir 3 )y)ri + 2r 2 (d 32 (w + 2r 3 (r 2 r 2 - {r\ +x- l)n + r 2 u) - 2(2r 2 + u)y) 
+d 21 (u + 2(r 3 r\ - r 2 r 3 r\ + r 3 (u + x- l)n + r\y - u(x + y))) + d 22 (u + 2(r 3 r\ 

-r 2 r 3 r\ + r 3 (u + x- l)n + r\y - u(x + y))))n + di(ri - r 2 )r 3 (r 2 - r 2 r x + u)), (B27) 

mz 1 N ° (_2d A ( ri - r 2 )r 3 (r 2 + r\ + u) - Sd A1 r ir2 {2r\ + r 2 r\ + Zur\ + r£n + u(u 
8Lirir 2y /ri 

-r 2 )) + 2d 21 r 1 r 2 r 3 (2r 1 x - r 3 ) + 2d 32 r 1 r 2 r 3 (r 3 - 2r 2 x) - d 1 (r 1 - r 2 )r 3 (2r 1 (r 1 - r 2 ) + 2x 
+2y - 1) + 4d 3 irir 2 (u + 2r 3 (r 2 r 2 - {r\ + x - l)n + r 2 u) - 2{2r\ + u)y) - d 2 r 3 {2r 3 r\ 
+(-2x -2y + l)n -r 2 + 2r 2 (x + y)) + d 3 (2r 3 y - 2r 3 yr x + r 2 (r 2 {2r\ - 1) + r\)) 
-Ad 22 rir 2 {2yrl + u + 2{r\ + (-r 3 + r 2 (u + x) + r\{—r\ + u + x))r ± - u(x + y))) 
-4diinr 2 (-« + 2(u - 2nr 3 )y + 2n(3r 3 r 2 - r 2 r 3 n + r 3 + (3n + r 2 )(u - x)) j), (B28) 
m z 7/2 N Jr 3 ~ 



(ri-r 2 )(diT 6 + dsT c ) (B29) 
r\r 2 L\ 



4. Spin-triplet P-wave state: |(c&)[ 3 Pj]) 

There are nine new Lorentz structures for |(c6)[ 3 Pj]), 
which are 
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B 1 = ^-e J aa k v , B 2 = -^k a k^ei^, B 3 = -K-k a <ge J afl k 1 ', B 4 = -Lq^e J a0 k\ 
tllZ z z z 

B b = —e(k, q 2 ,P, a)e J ap k v , B 6 = —e(q 2 , q 3 , P, a)e J afj k v , B 7 = —e{k, q 3 ,P, a)e J aj3 k v , 

777 ^ TTt ^ TTL ^ 

B 8 = -^ T e(k,q 2 ,q 3 ,a)k f3 ei l3 k u ,B g = -^ T e(k,q 2 ,q 3 ,a)q^ei p k' y , 

777/ ry 777 rj 



It is noted that e°'| is the symmetric tensor and e^g A*(|(cfc)i[ 3 -P 2 ])) =0 for i = (1 — 4), (B32) 

is the anti-symmetric tensor, and the fact that s x aa = j = (1,5,6,7,8,9). 

£ aa = 0. As the terms involving the following coefficients 
have no contributions to the square of the amplitude. So 
practically we can safely set the coefficients before them 
to be zero: 

4(|(c6)i[ 3 Po]» = for i = (1-4), 
j= (5,6,7,8,9) 

4(|(c6) 1 [ 3 P 1 ])) = fori = (l-4), mol , f The coefficients A], A 3 2 ,Af and that are the same 
3 (B31) for all three P- waves: 

3 = (1) 

I 



IT m ^ I 2 

A\ = 1 % (2d 21 r 2 T b (r 3 - 2r x x)r x 2 + 2d 32 r 2 T c (r 3 - 2r 2 x)n 2 + d 4 (3n + r 2 )r 3 T c 



x(rir 3 + u) + d 1 T b ((3r 1 + r 2 )r 3 (r 1 r 3 + u) - 2r 1 (r 1 - r 2 )x) - d 2 r 3 T b (r 1 (3r 1 2 + r 2 

+6r 2ri ) + (3n + r 2 )u - 2r 3 y) + d 3 r 3 T c (2n 3 + (r 3 2 + 3u)n + r 2 u - 2r 3 y)), (B33) 

Al = -16dnm z 5 / 2 ri x /^T 6 Li, (B34) 

A\ = -8mz 5 / 2 r 3 3/2 (rfnT 6 - d 22 T b + (d 31 + d 4X )T c )L u (B35) 

rir 2 

= WZ 7 fe (((3di + 3d 2 + cfe - d 4 )r 2 + 2(2di - d 2 + 2d 2 i + d 3 + 2d 32 )r 2 n 
2rir 2 r 3 ' 

+(di - rf 2 + 5d 3 + cLt)^)^ - 8rir 2 (d 2 iri + d 32 r 2 )x), (B37) 
Al = _ mz 5 / 2 ^i 

7-2 
\A3 



4= 4mZ !. 2Ll (fe-rf3i), (B40) 

2 777 Z^ ^ ^ 

-4? = » -=(2d 2 ir 2 (r 3 - 2r x x)ri 2 + 2d 32 r 2 (r 3 - 2r 2 x)r 1 2 + d 4 (3ri + r 2 )r 3 (rir 3 + u) 

2r Y r 2 Jri 



+d 1 (3r 1 4 + 7r 2 ri 3 + (5r 2 2 +u- 2x)r 1 2 + r 2 (r 2 2 - 8u + 2x)n - r 2 2 u) + d 2 (3r 3 ri 3 
+uri 2 + r 2 (6u - r 2 r 3 )ri + r 2 2 u - 2r 3 2 y) + d 3 (3r 3 ri 3 + 3wr x 2 + 2r 2 wi 

-r 2 2 r 3 n + 3r 2 2 u - 2r 3 2 y)), (B41) 
A\ = -U lx m z h l 2 r x ^- 3 L 2 , (B42) 
A\ = -2m z 5 / 2 ((3dn - d 22 + d 31 + d 41 )n - (d n - d 22 + d 31 - d 41 )r 2 )y^L 2 , (B43) 

A 2 = - Amz5 ' 2 ^ L \ d 2 T b + d 3 T c )( ri r 2 ), (B44) 
nr 2 



2 m 5/2r 

^-^(((di - d 2 )rj + 4(d 2 - d 21 )r 2ri - (di - d 2 )r 2 )r 3 T b + ((d 3 - d 4 )rf 

rir 2 r 3 / 

-4(d 3 - d 32 )r 2 ri - (d 3 - d 4 )r 2 )r 3 T c + 8rir 2 (d 2 iriT b - d 32 r 2 T c )x), 
Am z 5 / 2 J^T b L 2 



r 2 

16mz 5/2 L 2 



lr 3 

16m z 5 / 2 L 2 



'-(di - d 2 ), 
((d 2 i+d 22 )T b -d 32 T c ), 
(d 22 T b + d 31 T c ), 



1 -(m z 7/2 N (^-(r 1 (2rl + hr 2 r Y + r\)r\ + (3n + r 2 )wr| + (n - r 2 )u 2 ) 



8^r3~L 2 rir 2 
-2d 2 i(rir 3 + u)(2nx - r 3 ) + 2d 32 (rir 3 + u)(r 3 - 2r 2 x) + -^-(2(r 2 + x - l)rf 
+4r 2 (r| + x - l)r 2 + (2(r| + u + x - l)r| - u + 4r 2 y + 2uj/)n + r 2 (2ur| + u 
+2(5r| + 3u - 3x - 2y + l)r? + 2r 2 (r 2 ! - 4rfy - 2m/)) + — (6r? + 14r 2 r? 

nr 2 

+4u - 4x - 4y + 2)r 2 + (2(u - x - 2y + l)r 2 + u(2y - l))n + r 2 w(l - 2j/)) 
+-^ 2 -(2rf + 4r 2 rf + 2(r 2 + u + x + y - l)rf + 2r 2 (u + 2x + y- 2)r\ 

nr 2 

+(u + 2(r 2 (x - y - 1) - u(x + y)))n - r 2 (2yr\ + u- 2ux - 2uy)))), 

diim z 7/2 N ^/r 3 ~ . 

— {nr 3 +u), 



L 2 

mz V N °Vr3( d22 _ + dll ( 2rir3 + 2x + 2y - 1) + 2(d 4 i(r 2 + u) 



2L 2 

-d 22 (nr 3 -x-y) + d 31 (r 2 r 3 + y))) 
m z 7/2 N 0y ^ / , ^ 2 2 2 



rir 2 L 



2 



(-diT 6 r^ + d 2 (rf - + r 2 )T 6 + d 3 (n - r 2 )r 3 T c ), 



1 -= 1 -{m z 7/2 N a {Ad 21 r 2 T b xr 2 1 - 2d 2l r 2 r 3 T b r 1 - 2d 4 r2r c n - 2d 3 r 2 r 3 T c r x 



+2d 32 r 2 r 3 T c ri - 4d 32 r 2 T c xri - d 3 r 3 T c - 2d<±r 3 T c u + d 2 r 3 T b (2nr 2 — 2x — 2y + l) 
+2d 3 r 3 T c y + dir 3 T 6 (4x + 2y- 1))), 

_ m z 7 ' Np/ri ^ ri ^ _ r2)Tc _ 2d 2 T b (nr 2 + u) + d 1 T b (r 1 r 3 + 2u)), 
Am z 7 / 2 N Jn 



L 2 

±m z 7/2 N Jri 



(dnTb - d 32 T c ), 
{d u T b + d 31 T c ), 

—L-(m Z 7 / 2 N (- 2d * T ^ r * + H) (ri(2r § + U ) - , 2U ) + 2d 21 T b (r 2 - r 2 n + u) 
Zy^Li r x r 2 

x (2nx - r 3 ) - 2d 32 T c (r 2 - r 2 r Y + u)(r 3 - 2r 2 x) + — (2r 2 r? - 2(x - \)r\ 

r\r 2 

-2r 2 (r\ -u + 2x- 2)r\ + (2(u -x+ \)r\ + u - 2{2r\ + u)y)n + r 2 u{2y - 1)) 
'-(2(x + 2y - l)r\ + 4r 2 (3x + 2y - l)r 2 + (2(x + 2y - l)r 2 + u(2y - l))n 



nr 2 

d T 

+r 2 u(l - 2y)) + -^{2r 2 r\ + 2(x + y- l)r? + 2r 2 (-r 2 ! + u + 2x + 3y - 2)r\ 
r\r 2 " 

+(2(u + x + 3y - l)r| + u(2x + 2y - l))n + r 2 (2yr^ + it - 2mx - 2uy)))), 
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4 / 4d n m z 7 / 2 iVoV^r 6 2 
^2= ; W -r 2 ri (B58) 

- 2mZ? r 2Ar0x/ ^ ((rfnT b - rf 22 T b )(2x + 2y - 1) + T c (d 4 i(r 2 - r 2 2 + r 2 + 2«) 

+d 3 i(2y-l))), (B59) 

At = mZ l]^ [{d2 dl)rl {d2 d3)in r2)2) ' (B60) 

Af = 1 -—^{m z 7 ' 2 N {2d 21 r 1 {?,r 1 - r 2 )r 2 {2r lX - r 3 ) - 2d 32 n(r 1 - 3r 2 )r 2 (r 3 - 2r 2 x) 

8rir 2 r 3 ' L x 

+r 3 (dir 3 + 2{d 4 {r 2 1 + r\ + u) - c?i(r 2 - r 2 r x + 2x + y))r 3 + d 3 (2r 3 - 6r x r 2 

+r 3 - 2r 3 y) - rf 2 (ri(2n(n - 3r 2 ) + 1) + r 2 - 2r 3 x - 2r 3 y)))), (B61) 
= ^^^ ((rfi - 2^2 + ds)ri(n - r 2 ) + 2(di - d 2 )u), (B62) 
A|' = mz !!_ 2N ° ((rf 21 + d 22 + rf 32 )( ri - r 2 ) + (-du + rf 21 + rf 22 )r 3 ), (B63) 

4' = mZ !!_ 2 T N ° ((d 31 - d 22 )(n - r 2 ) + (du - d 22 )r 3 ), (B64) 

I 

Remaining non-zero coefficients for |(cfe)[ 3 Po]) an( l 
|(c6)PP 2 ]): 

I 



A\ = 8L 1 m z 5 / 2 x /^(rfii(3ri + r 2 )T b - r 3 ((rf 2 i + d 22 )T 6 + d 32 T c )), (B65) 
A 2 = 2L 2 m z 5 / 2 ((-di 1 + d 21 + d 22 + d 32 )r 2 - {-5d n + d 21 + d 22 + d 32 )n)^3, (B66) 
,,/ mz 7 / 2 NoTjr 3 

Af = — — - — (2rir 3 d 2 i - 2xd 21 + d 21 + d 22 + d 32 + 2d 22 r Y r 3 - 2d 32 r 2 r 3 - 2d 22 x 

IL 2 

-2(d 21 + d 22 + d 32 )y + du(4nr 3 + 2u + 2x + 2y - 1)), (B67) 

A A 3 = ~ 2mzV2 r Noy ^ (d 32 T c - (dai + d 22 )T b (2x - 1) - 2((d 21 + d 22 )T b + d 32 T c )y 

+d 11 T b (2r 1 (r 1 - r 2 ) + 2u + 2x + 2y - 1)), (B68) 

I 

Remaining non-zero coefficients for |(cfe)[ 3 Pi]): 



A\ = 8Lim z 5 / 2 V? 3 "(d 11 (r 1 - r 2 )T b + (d 21 + d 22 )r 3 T b + d 32 r 3 T c ), (B69) 
A\ = -2(dn - d 21 - rf 22 - d 32 )L 2 m z 5 l 2 ( ri - r 2 )V^, (B70) 

A\ = - mZ nT Noy ^ ' {r x r 2 {d 21 + d 22 + d 32 + 2((d 21 + d 22 )r 1 - d 32 r 2 )r 3 - 2(d 21 + d 22 )x 
AL 2 r\T 2 

-2(d 21 + d 22 + d 32 )y + dn(-2u + 2x + 2y- 1)) - d ir 2 ), (B71) 

A A 3 = - 2mZ l /2N ° s/¥i {d l T b r 2 + nr 2 (-d 32 T c + (d 21 + d 22 )T b (2x - 1) 

+d 11 T b {2r 1 (r 1 - r 2 ) + 2u - 2x - 2y + 1) + 2((rf 2 i + d 22 )T b + d 32 T c )y)), (B72) 

I 

Appendix C: Lorentz structure coefficients of L^ r , C 2 ~ r ' = ~ 7Uz ^ 1 ^ _|_ 4.J , c ^p 1 . g 3 (C2) 

All the non-zero coefficients of L^ r , , 

C ^ T ' = W (1 + 4T ^ 3 (C1) 
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CT r ' = ^^(1 + m)pi ■ P2 (C3) C 2 — ' = - 4(2 - V2)T d ) Pl ■ q 3 (C6) 

CT r,= X£f 1 (1 + 4Td) (C4) C 8 -^ = ^^(-V2 + 4(2-V2)T cJ )( Pl .p 2 ) (C7) 

Cr rr ' = ^^"(^2 - 4(2 - V2)T d )p 2 • g 3 (C5) CT^ = + 4T d - V2), (C8) 

where M = , 1 and Tw = i — sin 2 6L. Ni can be 

obtained by ensuring /cf = ^\^ >1V( " 7 V\vHz P V2a to satisfy 
/••-' 1. 
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